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ABSTRACT 


A criterion  of  fracture  under  proportional  straining  has  been  derived 
for  materials  that  undergo  plastic  deformation  prior  to  failure.  This  has 
been  done  using  endochronic  theories  of  plasticity  and  fracture  proposed  by 
Valanis.  The  strain  criterion  for  fracture  of  brittle  materials,  proposed 
by  Wu,  falls  into  the  class  of  more  specific  criteria  encompassed  by  the 
criterion  proposed  here. 

The  criterion  suggests  that  fracture  will  occur  when  the  magnitude  s of 

the  strain  vector  in  strain  space  reaches  a critical  value.  More  specifically, 

' 

in  the  event  that  s depends  only  on  the  angle  + between  the  strain  vector 

r ' 

and  hydrostatic  strain  vector,  it  then  follows  that  s ■ s($),  thus  giving 
rise  to  a polar  form  of  a fracture  surface  which  is  rotationally  symmetric 
with  respect  to  hydrostatic  strain  vector.  This  polar  form  allows  for  a one 
parameter  representation  of  the  fracture  surface  in  a triaxlal  strain  field. 

In  Part  II  of  this  paper  we  show  that  the  theory  is  fully  corroborated 
by  extensive  experimental  data  on  gray  cast  iron. 


□ □ 


I.  INTRODUCTION 


Fracture  of  materials  under  proportional  loading  has  been  a frequent 

/ 

topic  of  Investigation  since  it  plays  an  important  role  In  the  establishment 
of  safe  working  loads  In  engineering  design.  Numerous  fracture  criteria  have 
been  proposed  in  the  past.  A summary  of  the  available  theories,  though  not 
exhaustive,  can  be  found  in  Nadai  [1]  and  Paul  [2]. 

In  this  treatment,  we  have  taken  a new  approach  in  obtaining  a solution 
for  this  old  problem.  This  approach  stems  from  the  previous  works  of  the 
authors. 

An  energy-probability  theory  of  fracture  has  recently  been  proposed  by 
Valanis  [3].  This  theory  is  based  on  the  concepts  of  probability  and 
intrinsic  time;  the  latter  was  originally  introduced  by  the  first  author  in 
the  development  of  the  endochronic  theory  of  viscoplasticity  [4],  According 
to  the  aforementioned  theory,  fracture  occurs  on  an  intrinsic  time  scale  C 
which,  in  the  case  of  plasticity  is  given  by: 

d<2 ' Wk*"u  «•» 


where  1®  a positive  semi-def initc  rymmetric  material  tensor  which 

may  depend  on  the  Green  teneor  . Fracture  of  a microelement  will  have 
occurred  if  the  intrinsic  time  c of  the  microelement  has  reached  a critical 


value  cc»  which  is  given  by  the  equation: 


(1.2) 
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where  y is  a material  parameter;  k is  the  Boltzmann  constant;  T the 
absolute  temperature;  and  ^ is  the  change  in  energy  of  the  microelement 
(relative  to  its  unstressed  state).  The  fracture  activation  energy  ipQ 
introduced  in  Ref.  [3]  has  been  assumed  to  be  zero  in  the  present 
discussion. 

In  the  same  paper,  Valanis  discussed  the  case  of  fracture  under  pro- 
portional straining  when  the  material  obeys  Hooke's  law.  A brief  account 
of  this  discussion  is  given  in  the  following: 

In  the  Euclidean  principal  strain  space  (e^,  z 3)  let: 

ds2  - d£l  de1  (1.3) 

Proportional  straining  is  then  defined  by  the  relation 

dz±  - ds  (1.4) 

where  1 ± are  constant  direction  cosines  of  a radial  stain  path  of  length  s. 

Under  conditions  of  small  deformation,  the  Green  tensor  in  equation 

(1*1)  is  replaced  by  the  small  strain  tensor  . This  equation  in  con- 
junction with  equation  (1.4),  reduces  now  to 

dC2  - ^Vj  dfl2  (1.5) 

In  Ref.  (3]  a more  restrictive  assumption  was  made,  for  simplicity,  that 
depends  only  on  the  direction  of  the  vector  z±  but  not  on  its 
magnitude,  in  which  event,  the  scalar  product  1®  a function  of 

only.  In  addition,  a function  f(i  ) was  defined  such  that 


I 
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f - p l l 
ij  i J 


(1.6) 


Equation  (1.5)  may  therefore  be  written  as 


d<;  - ds 


(1.7> 


In  the  case  of  elastic  materials  obeying  Hooke's  law,  the  strain  energy, 
which  is  the  free  energy  at  constant  temperature,  is  given  by 


l|l  ■ 09 


(1.8) 


where  a " y and  cjj  *8  t*ie  “airlx  elastic  constants  defined 

by  the  linear  relation  a j ■ cij£j*  w^ere  are  the  principal  components 

of  the  stress  tensor. 

Substituting  equation  (1.8)  into  equation  (1.2)  and  Integrating,  the 
following  transcendental  equation  is  obtained  for  the  critical  values  s^ 
of  s : 


fai}  jBc  ■ V4^erfVkT  8c['  1 


The  solution  of  this  equation  is  of  the  fc’-m: 


8c  " 8cai) 


(1.9) 


(1.10) 


Equation  (1.10)  defines,  in  strain  space,  & fracture  surface  on,  or  within , 
which  all  fracture  initiation  events  will  have  occurred. 

Valanis  then  considered  the  case  where  the  fracture  surface  possesses 
rotational  symmetry  with  respect  to  the  hydrostatic  strain  vector  (Fig.  1). 
For  this  type  of  fracture  surface,  sc  is  a function  only  of  the  angle 
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between  the  strain  vector  and  the  hydrostatic  strain  vector,  where 

f?:  . 

cos  i - 'Y  ■ — (A^  + t~2  + *3)  (1.11) 

2 

and  i.  - tre  , i2  ■ tr(e  ) . 

The  function  sc  depends  on  through  a(A^)  and  f(A^).  It  was  shown, 

however,  that  any  asymmetry  of  the  fracture  surface  with  respect  to  the  null 
hydrostatic  strain  plane:  $ c -j,  must  be  due  to  the  function  f(A^). 

In  this  paper,  the  same  approach  as  described  above  with  the  fracture 
surface  having  rotational  symmetry  about  the  hydrostatic  axis  has  been 
followed  to  investigate  fracture  of  a class  of  materials  which  undergo 
plastic  deformation  prior  to  fracture.  In  particular,  a fracture  criterion 
for  this  class  of  materials  has  been  derived.  It  is  shown  in  Part  II  of 
this  paper,  that  in  the  case  of  gray  cast  iron,  the  form  of  this  criterion 
is  identical  to  the  strain  criterion  of  fracture  proposed  by  Wu  [5]. 

The  appropriate  constitutive  equation  used  in  the  derivation,  is  given 
in  the  Appendix.  This  equation  becomes  linear  in  s,  as  s Increases 
beyond  the  "elastic  range",  though  the  magnitude  of  the  strain  field  remains 
small  in  the  accepted  kinematic  sense.  The  physical  consequence  in  simple 
tension,  for  instance,  is  that  the  stress-strain  curve  becomes  asymptotically 
linear,  for  larger  "plastic"  strains. 

The  asymptotic  (linear)  form  of  the  constitutive  equation  plays  a signi- 
ficant role  in  the  prediction  of  fracture  of  plastic  materials.  The  obtained 
solution  is  a lower  bound. 

In  the  final  section,  we  discuss  at  length  the  material  tensor 
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which  appears  in  equation  (1.5).  This  tensor  is  pivotal  in  the  prediction 
of  fracture  after  the  material  has  experienced  a more  general  deformation 
history.  It  will  be  shown  that,  under  the  assumption  of  isotropy  of  the 
tensor  in  the  null  hydrostatic  strain  plane,  the  fracture  surface  is 

rotationally  symmetric  about  the  hydrostatic  strain  vector.  This  result 
will  be  justified  by  experimental  data  on  gray  cast  iron  in  Part  II  of  this 
paper. 


II.  The  Free  Energy 

For  plastic  materials  under  proportional  straining,  an  expression  for 
the  free  energy  function  ip  is  obtained  in  this  section,  which,  in  the 
next  section,  is  substituted  into  equation  (1.2)  to  obtain  the  fracture 
criterion. 

In  the  case  of  small  deformation  and  under  suitable  smoothness 
assumptions  regarding  the  functional  dependence  of  <p  on  g and  qa, 
one  can  express  the  free  energy  density  in  a quadratic  form  in  e and  qa, 
where  qa  (a  * 1,  2,  ....  n)  are  n internal  state  variables.  In  an  one 
dimensional  stress  field,  a direct  comparison  with  a mechanical  model  shows 
that  the  appropriate  form  is: 

. a 

where  are  the  spring  constants  and  and  q”  are,  respectively, 

the  strain  and  internal  state  variables  corresponding  to  the  one-dimensional 


field 
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In  a three-dimensional  situation,  isotropy  allows  a simple  decomposition 
of  iJj  into  the  form: 


* * *H  + *D 


where  iji^  and  i(/D  are  the  parts  of  free  energy  associated  with  volumetric 
deformation  and  shear  deformation  respectively,  as  can  be  readily  seen  in 
Appendix  II  of  Ref.  [4].  We  thus  have: 


where 


= ^tr  I(?  “ P“)2] 

a 

"H  " 7 XX  (e  - q“)2 


a a 1 . a 

P “ s - 3 « q 


a „ a 
q - tr  q 


e “ e - y 6 e 


e * tr  e 


The  constants  and  are  the  spring  constants  corresponding  to  shear 

and  volumetric  deformations,  respectively. 

Using  equations  (2.2  - 2.5)  and  following  Ref.  [4],  the  rate  equations 
of  the  internal  state  variables  are  obtained  in  the  principal  strain  space 
in  the  following  form: 


p 


a a dq 
1 + b -r3-  “ K_ 


K q + b - K e 

o dz  a 


(a  not  summed) 


(2.6) 


, a , a dpl  _ 

2paPl  +b2dT*  2pa  ei 


(a  not  summed) 


(2.7) 


where  bQ  and  b^  are  the  volumetric  and  the  deviatoric  parts  of  the 
viscosity  tensor  respectively  and  z is  the  intrinsic  time  scale. 
Integration  of  equations  (2.6)  and  (2.7)  gives: 


q“  “ e (z) 


-I. 


z -A  (z-z')  . 

a 3e  , 

e — , dz' 


3z 


(2.8) 


and 


h " ei(2)  '/ 

Jo 


z -A(z-z')  3e 

a i . . 

e JT dz’ 


(2.9) 


where,  for  constant  Poisson's  ratio*, 


K 2 p 

A --SL  --a- 


a .a  , a 
bo  b2 


(2.10) 


One  can  combine  equations  (2.8)  and  ('’.a)  to  obtain: 


z -X  ( z-z ')  3e- 
_ a i 


3z 


- dz' 


(2.11) 


We  now  observe  that  since  for  proportional  straining  eA  “ J^s  and 
e ■ /3  s coS(j> , equation  (2.11)  can  be  reduced  to: 


* This  is  a reasonable  assumption  as  has  been  demonstrated  in  Ref.'s  [6] 
and  [ 7 ] . 
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<x  „ -a 
qi  “ *iq 


(2.12) 


(2.13) 


uatlon  (2.12)  shows  that  for  proportional  straining,  the  Internal  variables 


where,  for  constant  Posson's  ratio,  the  heredity  functions  are  given  by: 


where 


and 


p(z)  = yQG(z)  , K(z)  = KqG(z) 


(2.15a) 

(2.15b) 

(2.15c) 
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The  Integral  in  equation  (2.14)  can  now  be  expressed  in  closed  form  by 
use  of  equation  (1.7)  and  the  following  relation  previously  Introduced  by 
Valanis  In  Ref.  [4], 


z - j log  (1  + 8?) 


(2.16) 


where  B is  a constant.  The  resulting  expression  Is: 


+ 2vii,2t)  • 5G4(i^f£)f1 ' (l  + Bfe>  "“] 


(2.17) 


where 


n 


a 


(2.18) 


Equation  (2.17)  represents  the  free  energy  for  plastic  materials  under 
proportional  straining.  For  large  s,  the  free  energy  i|i  has  the 
asymptotic  representation  t|/  , where 


*A  ■ i(3Koco,z» + 2v  ""20  • X)  °«(ip5f4): 


or 


*A"  2 


2 2 
3K  cos  * + 2p  sin  6 
o o T 


nri±jfa.]2 
JL  Bfn  J 


where 


n an 


a 


(2.19) 


(2.20) 


(2.21) 


It  will  be  shown  in  the  Appendix  that  for  large  s,  the  asymptotic 
representation  of  the  constitutive  equation  under  proportional  straining 
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has  the  form: 


2p 


+ /3  K cos<j 

/I  ' 


1 + Bfs 

(5f  n 


(2.22) 


where 


n “ 1 


/y  fa 

n 


(2.23) 


a a 


Direct  substitution  of  equation  (2.22)  in  equation  (2.19)  gives  rise 
to  the  following  expression  for  large  s: 


EG 

-s 


a n 


a a 


2 *o 


(2.24) 


where 


is  given  by  equation  (2.25): 


K4'30' 


(2.25) 


and  has  the  physical  significance  of  the  recoverable  energy  which  would  be 
obtained  if  the  material  unloaded  in  a linear  elastic  fashion  from  the 
stressed  to  the  unstressed  state. 

As  a result  of  equation  (2.24)  we  have  the  following  consequence: 


The  free  energy  of  a plastic  material  under  proportional 
loading  and  large  s,  is  a quadratic  function  of  the 
stress . 


I 
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In  other  words,  if/^  is  a quadratic  function  of  o^. 


Also  one  can  show  that- 


23  a 23 


(2.26) 


This  may  be  shown  as  follows: 


From  equation  (2.15c)  and  since  G are  all  non-negative  it  follows 

a 


that  G < 1 for  all  a.  In  that  event, 
a — 


Eg  2 t — \ u 

r iE-t 

a La.  or  n 


(2.27) 


However,  as  a result  of  the  Bessel  Inequality: 


(2.28) 


Hence,  equation  (2.26)  follows  directly.  As  a result  of  Inequality  (2.26) 
we  have  the  following  consequence: 


> t|i 

rA  — a 

For  Instance  in  simple  tension: 


(2.29) 


o 


(2.30) 


which  is  a revealing  result.  Note  that  for  one  Internal  variable  if/,  is 

A 


1 2 

precisely  equal  to  y o /Eq. 


! 

i 
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III.  The  Fracture  Criterion 

Substitution  of  equation  (2.17)  in  equation  (1.2)  yields  a fracture 
criterion  which  has  a representation  of  the  type: 

F(*,sc)  - 1 (3.1) 

from  which  s of  s at  fracture  may  be  determined.  This  fracture  criterion 

c 

has  a number  of  features  which  are  worth  discussing.  It  is  of  a polar  form 
which  can  be  written  as: 

sc  - F*(*)  (3.2) 


The  significance  of  equation  (3.2)  is  that  it  is  an  one  parameter  represen- 
tation of  the  fracture  surface  in  the  presence  of  a three  dimensional  strain 
field.  It  is  analytically  defined  in  terms  of  a small  number  of  parameters, 
and  a function  f (<(>). 

In  the  particular  case  where  fracture  occurs  for  a sufficiently  large 

s (in  simple  tension  this  would  be  on  the  straight  part  of  the  "plastic" 

portion  of  the  stress-strain  curve) , then  the  asymptotic  expression  ^ of 

A 

given  by  equation  (2.20)  may  be  used  as  an  approximation;  this  leads  to 
a closed  form  solution.  Substitution  of  equation  (2.20)  in  equation  (1.2), 
then  gives  the  following  expression  for  8c. 


(3.3) 
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where  v is  Poisson's  ratio. 

It  may  be  noted  that  the  free  energy  is,  thus,  continuously  overestimated 
and  the  critical  value  sc  determined  from  the  resulting  equation  (3.3)  is, 
in  fact  a LOWER  BOUND. 

In  Part  II  of  this  paper  we  show  that  the  lower  bound  is  a good 
approximation  to  the  actual  value  as  far  as  gray  cast  iron  is  concerned. 

Equation  (3.3)  has  been  derived  on  the  basis  of  the  assumption  that  the 
fracture  surface  is  rotationally  symmetric  about  the  hydrostatic  strain  axis. 
This  assumption  will  be  further  discussed  later.  We  also  notice  from  equation 
(3.3)  that  any  asymmetry  of  the  fracture  surface  with  respect  to  the  null 
hydrostatic  strain  plane  ij>  ■ y must  be  due  to  the  function  f(i^).  This 
conclusion  has  been  obtained  for  elastic  materials  by  Valanis  [3]  and  has 
been  shown  here  to  be  valid  for  plastic  materials  as  well. 

Symmetries  of  the  Fracture  Surface 

The  fracture  locus  in  the  deviatoric  plane  (<f>  = y ) must  possess  certain 
symmetries  for  isotropic  materials  on  the  basis  of  the  following  considerations 

(1)  With  veference  to  Figure  1,  the  shape  of  the  fracture  locus  should 
remain  invariant  with  permutation  of  the  indices  1,  2 and  3 of  the  strain-axes. 
The  cons:raint  resulting  from  this  consideration  is  that  the  fracture  locus 
must  possess  a three-fold  symmetry  with  respect  to  the  axes  e^,  and  e^. 

(2)  It  is  observed  that  the  fracture  properties  remain  invariant  with 
a reversal  of  the  sign  of  the  shear  field  in  the  presence  of  uniform  hydro- 
static strain.  The  implication  of  this  consideration  is  that  the  fracture 
locus  in  the  deviatoric  strain  plane  should  be  symmetric  with  respect  to  the 
point  of  intersection  of  this  plane  and  the  hydrostatic  axis. 
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The  first  consideration  requires  that,  in  so  far  as  isotropic  materials 

are  concerned,  a proposed  fracture  surface  must  be  symmetric  with  respect  to 

the  axes  e^,  e2  and  the  second  consideration  requires  that  it  is  also 

symmetric  with  respect  to  the  axes  e j , and  which  are  orthogonal  to 

the  axes  e^,  e 2 and  e3  respectively,  at  0.  Therefore,  OAj  ■ OA^, 

0A2  » QA2  and  OA^  - 0A^  in  Figure  1.  Combination  of  the  above  two  consi- 
* 

derations  requires  of  the  fracture  surface  a six-fold  symmetry  with  respect 
to  the  six  axes  e^,  c2,  e^,  e2»  £3-  Thus,  only  a 30°  sector  of  the 

fracture  locus  need  to  be  found  for  the  entire  locus  to  be  completely 
determined . 


Material  Tensor  P, 


If  the  matrix  P were  known,  f would  be  completely  determined  in 
consequence  of  the  relation  (1.6),  and  the  fracture  surface  would,  therefore, 
also  be  determined.  Conversely,  equation  (3.1)  contributes  toward  the  deter- 
mination of  P^j . 

The  material  matrix  P^  is  now  shown  to  be  a second  order  tensor  in  the 

tensor  space  t^.  This  follows  from  the  axiom  below: 

2 

Axiom:  The  functional  dependence  of  d£  on  de^  is  form- invariant 

with  respect  to  rotation  of  the  coordinate  axes  in  the  tensor 
space  ei>  i.e.. 


d s2  - dEi  dej  ' dei  dEj 


(3.4) 


These  arguments  are  reminiscent  of  similar  ones  in  conjunction  with  the 
yield  locus  in  the  classical  theory  of  plasticity.  There,  incompressibility 
and  absence  of  Bauschinger  effect  were  assumed.  Such  considerations  do  not 
arise  in  this  discussion. 
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where 


d£i  - Qy  dej 


0.5) 


and  Is  a proper  orthogonal  transformation. 

This  axiom  is  then  a necessary  and  sufficient  condition  for  to  be  a 

tensor.  To  prove  the  sufficient  condition,  we  note  that  equations  (3.4)  and 
(3.5)  lead  to  the  expression 


<Fk«  - p;AkV’  d€k  dt«  - ° 


(3.6) 


Since  de^  is  an  arbitrary  vector  and  P is  symmetric  it  follows  that 


Pki  “ QikQjlPij 


(3.6) 


which  establishes  the  aforementioned  tensorial  property  of  P . 

Conversely,  if  P is  a second  order  tensor,  then  substitution  of  (3.5) 
and  (3.7)  into  (3.4)  gives 


dc  - Ptj  dEi  dEj  - QklQAJP^Qlr  de;  qj8  dc; 


6.  6 P'  de ' de'  - p;  „ tie'  de! 
kr  is  k£  r s ki  k i 


(3.8) 


Thus,  the  above  axiom  is  recovered. 

We  remark  that  the  above  proof  that  P^  is  a tensor  in  the  tensor 
space  eA  follows  directly  from  the  definition  of  the  endochronlc  time  given 
by  equation  (1.1),  and  does  not  involve  the  fracture  properties  of  a material. 
Furthermore,  the  above  proof  may  be  easily  generalized  to  include  the  case 


of  six-dimensional  tensor  space  e 


It  is  convenient  to  introduce  a primed  orthogonal  coordinate  system  such 
that  3'-axis  is  also  the  axis  of  the  hydrostatic  strain  (see  Figure  2). 
Symmetries  of  can  now  be  deduced  from  the  6-fold  symmetry  of  the  fracture 

locus  in  the  deviatoric  plane.  The  latter  also  implies  a 6-fold  symmetry  for 
the  function  f(Jl^)  on  the  same  plane.  Thus,  the  value  of  f i6  the  same 
at  points  A,  B,  C,  D,  E and  F.  In  the  particular  case,  when  the  fracture 
surface  is  rotationally  symmetric  with  respect  to  the  hydrostatic  axis,  the 
value  of  f remains  the  same  everywhere  on  the  circle  ABCDEF.  Let  be 

the  components  of  tensor  P referred  to  the  primed  axes;  let  and  be 

the  direction  cosines  of  two  proportional  strain  paths  that  correspond  to  the 


same  value  of  f(i|).  These  direction  cosines  are  related  by 


*i  - Vi 


(3.9) 


where 


cos0  sin0  0 


= -sin0  cos0  0 


(3.10) 


and  0 is  the  angle  between  the  projections  of  the  vectors  and  in  the 

deviatoric  strain  plane.  For  the  case  of  6-fold  symmetry,  0 is  a multiple 

* if 
of  -j. 

2 

For  the  class  of  vectors  l"  for  which  f remains  constant  we  have 


the  relation 


f z m pt  0 * 0 I . pi  ,11.11 

ij  i j ij*i  j 


(3.11) 
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Substitution  of  equation  (3.9)  In  equation  (3.11)  yields  readily  the  important 
relation 


P!  4 - R4,  R^P^, 
ij  ik  J1  ki 


(3.12) 


which  defines  the  symmetries  of  the  tensor  P. 

In  particular,  substitution  of  equation  (3.10)  in  equation  (3.12)  and  in 
the  event  that  0 Is  a multiple  of  y,  it  transpires  that,  for  these  values 
of  9 : 


p*  m p'  4 p' 

11  22  T 33  ’ 


p ' * p'  « p ' « 0 

r12  r23  31 


(3.13) 


Therefore,  at  the  axes  of  symmetry,  the  tensor  P^  has  only  two  independent 
components.  The  off-diagonal  terms  are  zero.  This  does  not  imply,  however, 
that  p^2’  P23  an<*  P13  are  zero  at  ot^er  locations  in  the  deviatoric  strain 
plane.  In  fact,  all  components  of  P^  are  periodic  in  magnitude  with  an 
angle  of  periodicity  y.  Furthermore,  P^  is  positive  semi  definite.  P^  also 
enjoys  the  reflective  symmetry  with  respect  to  the  axes  e^,  and 
(see  Figure  1).  This  implies  that  in  determining  the  fracture  locus  of  an 
isotropic  material  experimentally,  only  a 30  degree  sector  need  be 
Investigated  as  stated  previously. 

In  the  particular  case  where  the  fracture  locus  in  the  deviatoric  plane 
is  circular  then  equation  (3.12)  is  true  for  all  defined  by  equation 

(3.10),  i.e.,  for  all  0.  In  this  event  P^  is  Independent  of  0 and 

f(*i)2  " “ p[i<42  + ^2)  + ?33l32  “ Pii8in2*  + P33cos2*  “ 

(3.14) 


m 


Ld 


18 


That  is,  f is  only  a function  of  $ . If  is  fixed,  then  f is  a 

constant.  Conversely,  if  P^  is  independent  of  0,  then  the  fracture  locus 
in  the  deviatoric  strain  plane  is  a circle. 

Finally,  we  show  that  the  material  tensor  P can  not  be  fully  isotropic. 
This  can  be  shown  by  assuming  P^  to  be  isotropic  and  then  show  that  this 
leads  to  a consequence  which  is  not  in  agreement  with  observation.  For  isotropic 
P^j , we  can  write 


Pij(ek) " * V8’  W 


The  most  general  isotropic  of  P . form  is 

**  J 


Pij  “ Po(s)6ij  + Pl(s)Lij  + P2(s)LikLkj 


where 


, tr  L » tr  = tr  « 1 


Equation  (3.16)  can  further  be  simplified  to  obtain 


Pij  “ Po6ij  + (Pl+P2)*’i*’J 


Substitution  of  (3.18)  in  (1.6)  yields 


f -Vi*j  ”Po  + Pl  + P2 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


1 


« a function  of  s only 


(3.19) 


Thus,  f is  Independent  of  $.  In  particular,  if  P^  is  a function 
of  only,  then  f is  also  independent  of  s.  Therefore,  as  the  result 


m 


of  fracture  criterion  (3.3),  the  fracture  surface  must  be  symmetric  with 
respect  to  the  deviatoric  plane  in  contradiction  with  observation. 

It  is  interesting  to  remark  that  the  assumption  of  transverse  isotropy 
of  with  respect  to  the  hydrostatic  strain  axis  leads  readily  to  the 

following  result,  though  the  application  of  equation  (3.6): 


1-M  M 


Pjj  “ N M 1-M  M 

M M 1-M 


(3.20) 


where 


M - - 


(pi  - ri) 


(P[  + 2P^  * N " 3 (P1  + 2P3)  * 


(3.21) 


In  the  above  equations,  Pj  and  P^  are  the  two  independent  parameters  in  Pj  , 
We  note  that  P is  thus  determined  by  two  independent  parameters  M and  N. 

This  last  form  (3.20)  also  follows  if  P ^ is  of  the  special 
isotropic  form: 


Pijk*  “ Pl6iJ6kl  + 2 P2  6ik6U  + 6ii6jk 


(3.22) 


where,  however,  P and  P are  functions  of  cco  <j>  which  is  given 

by  eq.  (1.11);  in  other  words  P^  and  are  not  constants  but  are  functions 

of  a certain  sort  of  the  Invariants  of  e^j* 


IV.  APPENDIX  - THE  STRESS-STRAIN  RELATION  UNDER  PROPORTIONAL  STRAINING 


In  the  discussion  of  fracture  of  plastic  materials,  a stress-strain 


relation  is  useful.  This  can  be  derived  from  the  endochronic  theory  of 
plasticity,  in  which 
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„ _J±_ 


ij  3e 


(A.l) 


lj 


This  equation,  in  conjuction  with  equations  (2.2  - 2.4),  gives  rise  to  the 
following  equation  in  the  principal  space: 


«i  • - <$> + «.  - ! %>"  - <“>] 


which,  by  use  of  equations  (2.8)  and  (2.11),  can  be  reduced  to 


(A. 2) 


’i  - £{  - ^) + ^ v-*] /o 


z -X  (z-z' ) j 

e |f,  dz*  (A.3) 


or,  for  proportional  straining: 


°i  - [^(‘l  - ^f)+  * V0'*]  + ”“]| 

(A. 4) 

The  above  equation  governs  the  relation  between  the  stress  components  o ^ and 
the  strain  magnitude  s for  a given  angle  <f>.  For  large  s,  the  stress- 
strain  relation  (A. 4)  may  be  approximated  by  the  following  asymptotic  form: 


0 ^ “ C°j^  + *^3  K^cosiji 


/3  / 


A+_gfs\ 

V efn  ) 


(A. 5) 


where 


n - 1 


/z 


(A. 6) 


The  parameters  n and  8f  in  equation  (A. 5)  can  be  determined  from  the 
uniaxial  tension  test.  In  this  test,  one  has  a proportional  stress  path;  one 


i 
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also  has  a proportional  strain  path  if  Poisson's  ratio  is  a constant.  In 
particular. 


-v  e. 


(A.  7) 


1 

A + 2V2 


and 


Furthermore,  equation  (A. 4)  becomes 


(A. 8) 


(A.  9) 
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I 

I I 

where  E(o)  denotes  the  initial  slope  of  the  uniaxial  stress-strain  curve, 

is  the  tangent  modulus  at  large  strain  (the  slope  of  the  asymptotic 

straight  line),  and  a the  intercept  of  the  asymptotic  straight  line  with 

o 

the  stress  axis.  Equations  (A. 12)  and  (A. 13)  are  used  to  determine  the 
parameters  n and  for  a given  material. 
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FIGURE  CAPTIONS 

Figure  1 Fracture  locus  in  the  deviatoric  strain  plane  (<J>  ■ 
Figure  2 The  transformation  R^ 
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